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quasi-linear closure and an eddy-damped fluid model. Results indicate that despite 
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1. INTRODUCTION . 
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, , , . , . . , , sheet formation and, thus, reconnection (see 

Magnetic reconnection is the process 
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Alfven speed associated with B; S = vaL/t] 
is the Lundquist number; R m = ul/r) is 
the magnetic Reynolds number, with u 
and / being the characteristic amplitude 
and length scale of the velocity — S is 
called the magnetic Reynolds number 
R m in some literatures. Note that the 
SP process forms strangely anisotropic 
current sheets since A/L = y/S and 
S >> 1. Note also the link between sheet 
anisotropy and the reconnection speed v r , 
i.e., v t /va = A/L = 1/VS. Finally, it should 
be noted that v r is a measure of the global 
reconnection rate, in that it parameterizes 
the mean inflow velocity to the layer. 

The SP Picture is intrinsically appealing, 
on account of its simplicity and dependence 
only upon conservation laws. Moreover, the 
SP prediction has been verified by laboratory 
experiments (Ji, Yamada, & Kulsrud 1998). 
However, since R m is extremely large in 
most astrophysical applications of interest 
(i.e. R m ~ 10 13 in the solar corona), the 
SP reconnection speed is pathetically slow. 
Hence, there have been many attempts 
to develop models of fast reconnection. 



For example, in 1964 Petschek proposed 
a fast reconnection model involving shock 
formation near the reconnection layer, which 
predicted v r = va/ In S. Unfortunately, 
subsequent numerical (Biskamp 1986) 
and theoretical (Kulsrud 2000) study has 
indicated that Petschek's model is internally 
inconsistent. While research on fast, laminar 
reconnection continues today (i.e., Kleva, 
Drake, & Waelbroeck 1995) in the context of 
two-fluid models, the failure of the Petschek 
scenario has sparked increased interest 
in turbulent reconnection (Matthaeus & 
Lamkin 1986) in which turbulent transport 
coefficients (which can be large for large 
Reynolds number) act as effective dissipation 
coefficients, and so are thought to facilitate 
fast reconnection (i.e., Diamond et al. 
1984; Strauss 1988). Interest in turbulent 
reconnection has also been stimulated by 
the fact that many instances of reconnection 
occur in systems where turbulence is 
ubiquitous, i.e., coronal heating of turbulent 
accretion disks, the dynamo in the sun's 
convection zone, and turbulent tokamak 
plasmas during disruptions. 
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Recently, Lazarian and Vishniac (1999) 
(referred to hereafter as LV) presented a 
detailed discussion of turbulent reconnection. 
LV took a rather novel approach to the 
problem by considering the interaction of 
two slabs of oppositely directed, chaotic 
magnetic fields when advected together. 
LV modeled the effects of turbulence by 
treating the slabs' surfaces as rough, where 
the roughness was symptomatic of a chaotic 
turbulent magnetic field structure. This 
'rough surface' model naturally led LV to 
decompose the reconnection process into 
an ensemble of local, 'micro'-reconnection 
events, which interact to form a net 'global' 
reconnection process. LV argue that micro- 
reconnection events occur in small scale 
'layers', with dimensions set by the structure 
of the underlying Alfvenic MHD turbulence 
(i.e., the kj 1 and k^ 1 , as set by the 
Goldreich-Sridhar model). The upper bound 
for the micro-reconnection rate obtained 
by LV is v r = v A (u/v A ) 2 = v A (b/B H ) 2 , 
where Bh is the mean, reconnection field, 
and u and b are small-scale velocity and 
magnetic field. While the LV arguments 



concerning micro-reconnection are at least 
plausible, their assertion that the global 
reconnection rate can be obtained by 
effectively superposing micro-reconnection 
events is unsubstantiated and rather 
dubious, in that it neglects dynamical 
interactions between micro-layers. Such 
interactions are particularly important 
for enforcing topological conservation 
laws. Since the process of turbulent 
reconnection is intimately related to the 
rate of flux dissipation, and the latter 
is severely constrained by mean square 
magnetic potential conservation, it stands to 
reason that such a topological conservation 
law will also constrain the rate of global 
magnetic reconnection. In particular, for 
a mean _B-field with strength in excess of 
B ciit ~ y/(u 2 )/R m , the flux 2D was shown 
to be suppressed by a factor 



1 + R m (B}y (u 2 ) ' v 
where (B) is the large-scale magnetic 
field and (u 2 ) the turbulent kinetic energy 
(Cattaneo & Vainshtein 1991; Gruzinov 
& Diamond 1994); The above expression 
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implies that even a weak magnetic field 
(i.e, one far below the equipartition value 
(b 2 ) ~ (u 2 )) is potentially important. The 
origin of this suppression is ultimately 
linked to the conservation of mean square 
potential (see Das & Diamond 2000 for flux 
diffusion in EMHD). Hence, it is natural to 
investigate the effect of such constraints on 
reconnection, as well. 

In turbulent reconnection, fluctuating 
magnetic fields are dynamically coupled to a 
large-scale magnetic field so that a similar 
suppression of energy transfer is expected to 
occur. In other words, fluctuating magnetic 
fields will inhibit the energy transfer from 
large-scale to small-scale magnetic fields 
(responsible for turbulent diffusion), even 
when the latter is far below equipartition 
value. This link between small and large 
scale magnetic field dynamics is indeed the 
very feature that is missing in LV, where 
a global reconnection rate is considered 
to be a simple sum of local reconnection 
events, without depending on either (B) or 
R m . That is, even if one local reconnection 
event may proceeds fast, the energy transfer 



from large-scale to small-scale is suppressed 
inversely with R m , preventing many local 
reconnection events for a large R m and fixed 
large-scale field strength. Thus, the global 
reconnection rate is very likely to be reduced 
for large R m . 

The purpose of this paper is to 
determine the global reconnection rate 
by treating the dynamics of large and 
small-scale magnetic fields in a consistent 
way. The key idea is to compute the effective 
dissipation rate of a large-scale magnetic 
field (turbulent diffusivity) by taking into 
account small-scale field backreaction and 
then to use Sweek-Parker type balance 
relations to obtain the global reconnection 
rate. Since magnetic fields across current 
sheets are not always strictly antiparallel 
in real systems, we assume that only one 
component of the magnetic field (e.g., 
poloidal or horizontal field) changes its 
sign across the current sheet (see Fig. 2). 
The other component (e.g., axial field) is 
assumed to be very strong compared to 
the poloidal component. A strong axial 
magnetic field avoids the null point problem 
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inherent in SP slab model, justifying the 
assumption of incompressibility of the 
flow in the poloidal (horizontal) plane. 
Such a magnetic configuration is ideal for 
the application of so-called 3D reduced 
MHD (3D RMHD) (Strauss 1976). In 
3D RMHD, the conservation of the mean 
square potential is linearly broken due to 
the propagation of Alfven waves along an 
axial field, but preserved by the nonlinearity. 
As we shall show later, the latter effect 
introduces additional suppression in the 
effective dissipation of a large-scale magnetic 
field compared to 2D MHD. We also discuss 
the 2D MHD case which can be recovered 
from our results simply by taking the limit 
Bq — > 0, where B is a axial magnetic field. 

To be able to obtain analytic results, 
we adopt the following two methods. The 
first is a quasi-linear closure together using 
r approximation by assuming the same 
correlation time for fluctuating velocity and 
magnetic fields employing unity magnetic 
Prandtl number. The second is an eddy- 
damped fluid model, based on large viscosity 
(Kim 1999), which may have relevance in 



Galaxy where v ^> rj. In this model, the 
nonlinear backreaction can be incorporated 
consistently, without having to invoking the 
presence of fully developed MHD turbulence, 
or assumptions such as a quasi-linear closure 
or r approximation. In both models, the 
isotropy and homogeneity of turbulence is 
assumed in the horizontal (poloidal) plane 
since the reduction in effective dissipation 
of a large-scale poloidal magnetic field 
is likely to occur when its strength is far 
below the equipartition value. The effect 
of hyper-resistivity is incorporated in our 
analysis. This can potentially accelerate the 
dissipation of a large-scale poloidal magnetic 
field. 

The paper is organized in the following 
way. In §2, we set up our problem in 
3D RMHD and provide the quasi-linear 
closure using r approximation where the 
flux is estimated in a stationary case. 
Section 3 contains a similar analysis for 
an eddy-damped fluid model. The global 
reconnection rate for both models is 
presented in §4. Our main conclusion and 
discussion is found in §5. 
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2. QUASI-LINEAR MEAN FIELD 
EQUATIONS 

We assume that a strong constant 
axial magnetic field B is aligned in the z 
direction and that a poloidal (horizontal) 
magnetic field B# lies in the horizontal x-y 
plane, as shown in Fig. 2. The subscript 
H denotes horizontal direction. The 
total magnetic field is then expressed as 
B = Bqz + B// = B z + Vxt^i, in terms of 
a parallel component of the vector potential 
ip (i.e., B// = Vxipz). According to the 
RMHD ordering, the flow in the horizontal 
plane u is incompressible and therefore can 
be written using a scalar potential as 
u = Vx(f)z. Then, the equations governing 
3D RMHD are (see Strauss 1976): 

dt^ + u-Vip = r]V 2 ^ + B d z ^, (2) 
<9 t V 2 + u • VV 2 = z/V 2 V 2 + B • VV 2 ^3) 

where f] and v are Ohmic diffusivity and 
viscosity, respectively. For the quasi-linear 
closure, unity magnetic Prandtl number 
(77 = v) will implicitly be assumed. In 
comparison with 2D MHD, the equation for 
the vector potential contains an additional 



term B d z <f), which reflects the propagation 
of Alfven wave along the axial magnetic 
field B z. Due to this additional term, 
the conservation of the mean square 
potential is broken in 3D RMHD, albeit 
only linearly. In other words, the nonlinear 
term in equation (2) conserves (ip 2 ) since 
(u ■ Vip 2 ) = V • (uip 2 ) = 0, assuming that 
boundary terms vanish (cf. Blackman & 
Field 2000). Similarly, the momentum 
equation contains an additional term 
-Bo<9 2 V 2 t/>. These additional terms are 
proportional to the wavenumber k z along 
B z. Thus, the 2D case can be recovered 
by taking k z — > or B — > 0. Note that 
due to a strong axial field B z, the vertical 
wavenumber k z is much smaller than 
horizontal wavenumber k# = k x x + k y z; 
specifically, the 3D RMHD ordering implies 
that k z /k H ~ B H /B ~ e < 1. 

We envisage a situation where large- 
scale magnetic fields with a horizontal 
component B# = (B#) = Vx(ip)z are 
embedded in a turbulent background. 
The turbulence can be generated by an 
external forcing, for instance. The horizontal 



- 8- 



component of a large-scale magnetic field 
(Br) flows to form a current sheet of 
thickness A in the horizontal plane, so (B#) 
changes sign across the current sheet. As 
reconnection proceeds, small-scale flows as 
well as magnetic fields are generated within 
the current sheet. It is reasonable to model 
the physical processes within a current sheet 
as well as the background turbulence by an 
(approximately) isotropic and homogeneous 
turbulence with fluctuating velocity u 
and magnetic field b = Vx^'i Here the 
assumption of isotropy is justified since 
(B H ) 2 <C (u 2 ), i.e. the reconnecting field is 
taken to be weak. 

Outside the reconnection region, there 
are large-scale inflow and outflow in addition 
to the background turbulence. Thus, to 
obtain SP-like balance relations, small-scale 
flow as well as large-scale flow should be 
incorporated. However, since small-scale 
velocity is assumed to be homogeneous and 
isotropic, there is no net contribution from 
the fluctuating velocity to mass continuity. 
Effectively, the small-scale velocity does not 
appear in the momentum balance either. 



However, Ohm's law (magnetic energy 
balance) now contains an additional term 
due to the correlation between fluctuating 
fields (u x b), leading to turbulent diffusitivy 
(effective dissipation rate), which then 
effectively changes the Ohmic diffusivity to 
the sum of Ohmic diffusivity and turbulent 
diffusivity inside current sheet. Therefore, 
similar balance relations to the original 
SP hold in our case as long as the Ohmic 
diffusivity is replaced by the total diffusivity. 

To recapitulate, homogeneous and 
isotropic turbulence is assumed to be present 
with magnetic fields B# = (B H ) + h 
((b) = 0) and small-scale velocity u 
((u) = ((f)) = 0). Once the effective 
dissipation rate of (Bh) within the 
reconnection zone is computed, it will be 
used to determine the reconnection velocity 
v r through SP balance relations by using the 
total diffusivity in place of Ohmic diffusivity. 

2.1. Mean Field Equation 

The evolution equation for ip is obtained 
by taking the average of the above equation 
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as: 



following form. 



d t (iP) + (u-Vi/j') = W 2 ^). (4) 

Note that although equation (4) does not 
exhibit an explicit dependence on B , it does 
depend on B through the flux Tj = (uiip 1 ). 
To compute the flux T i: we first do a 
quasi-linear closure of (u ■ Vip'). 

The effect of the backreaction can be 
incorporated in the flux Tj by considering 
the change in flux Tj to be due to the change 
in the velocity as well as the fluctuating 
magnetic field. That is, we can rewrite the 
flux as 

r, = e^U') = Ujz&W - 8 08^%) 

where unity magnetic Prandtl number is 
assumed for the equal splitting between 
(dj4>8ip') and (S^djtp'); the latter essentially 
takes the backreaction to be as important as 
the kinematic contribution. 



2.2. Fluctuations 

^From equations (2) and (3), we can 
write the equation for the fluctuations in the 



(d t + u • V)^' - (u • VV0 = -u • V(^) + 77VV + L 
(d t + u- V) V 2 - (u • VV 2 0> = z/V 2 V 2 + B d z V V + 

Here we have assumed that there is no 
large-scale flow in the current sheet. To 
estimate 5(f) and 81])' in equation (5), 
we introduce a correlation time r that 
represents the overall effect of inertial and 
advection terms on the left hand side of the 
above equations. That is, we approximate 
(dt + u • V)V>' - (u • W) = r"V, and 
(dt + u • V)V 2 - (u • VV 2 0) = r^VV, 
where the same correlation time r is assumed 
for both the fluctuating flow and magnetic 
field due to unity magnetic Prandtl number. 
Then, 5(f) and 5ip' in equation (5) can be 
estimated from the above equations as 
follows: 

8i>' = T[B d z 4)' -e^d^'d^)] , 
<5V 2 = r[B d z V 2 ^' + e lj3 d 3 (^)d l V 2 ^' + e lj3 d^'d l V 2 (i 

In Fourier space, the above equations take 
the following form: 



'(k) = riBotk^ + e^Jd'k'k'^ik'^k-k^ik 



5i) 



5m 



IT 



BqWQs) + ie ij3 - / d 3 k' [(k - k') 3 k[k' 2 + 
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x^Ji^iQiift'^Jnnot be everywhere isotypic. 
Specifically, the correlation functions at 
equal time t are taken to have the form: 



Note that in principle, the correlation time 
can be a function of the spatial scale, or the 
wavenumber, i.e., r = t^. Nevertheless, for 
the notational simplicity, we have taken r to 
be a constant by assuming that the variation 
of Tk in k is small or that the small-scale 
fields possess a characteristic scale with 
a small spread in k. Our final result will 
not fundamentally change when the scale 
dependence of r is incorporated. 

The flux Ti can readily be computed 
once the statistics of small-scale magnetic 
field and the velocity are specified. As 
mentioned earlier, the statistics of both 
fluctuations are assumed to be homogeneous 
and isotropic in the x-y plane. We further 
assume that the former is homogeneous and 
reflectionally symmetric in the z direction 
with no cross correlation between horizontal 
and vertical components, thereby eliminating 
a helicity term. The absence of helicity terms 
rules out a possibility of a mean field dynamo 
in our model. Note that due to the presence 
of a strong axial field B z, the correlation 



(^(ki^'Ck^t)) = ^ki-kaMfciH.fcUp) 

((Kk!,*)^, t)> = 5(k 1 -k 2 )0(k 1H ,ka%) 

where tp(kiH, k lz ) and (p(k 1H , k lz ) are the 
power spectra of ip' and 0, respectively. 
These depend on only the magnitude of 
horizontal wavenumber kin = \J k\ x + k\ y 
and vertical wavenumber k lz . Finally, we 
assume that (00') = 0, which can be shown 
to be equivalent to excluding the generation 
of a large-scale flow by the Lorentz force. 

Straightforward but tedious algebra 
using equations (8)— (1 1) in equation (5) 
leads to the following expression for the flux 
(the details are given in Appendix A): 

r, = -\ W) - <& 2 ))W> - <V /2 >^v 2 <#]2) 

where (u 2 ) = j d 3 kk 2 0{k), (^ /2 ) = 
jd 3 k^(k), and (b 2 ) = jd 3 kk 2 ^{k). The 
first term on the right hand side of equation 
(12) represents the kinematic turbulent 
diffusion by fluid advection of the flux; 
the second represents the flux coalescence 
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due to the backreaction of small-scale 2.3. Stationary Case: d t (ib' 2 ) = 

magnetic fields with the (negative) diffusion 

rc ■ i , , n ii i To compute the flux I\, we need an 

coefficient proportional to the small-scale 

/l2\ rpi , i • i , ,i additional relation between (b 2 ) and r\ 

magnetic energy (o ). lhe third term is the x ' 

, . ,. ., a ,. , ,. besides equation (13). This can be attained 

hyper-resistivity, reflecting the contribution 

, t-, , i . r i i by imposing a stationarity condition on 

to f i due to the gradient ol a large-scale 

, T , _,. .. , T . _ _ N i^' 2 )- The stationarity of fluctuations is 

current (J) = -V 2 (ib). (Jz = VxB H ). XY 1 J 

it , ,i , i , c i . ,. ., achieved in a situation where the energy 

Note that the value of hyper-resistivity, 

, . , , i ■ i transfer from large-scale fields balances 

being proportional to mean square potential, 

„ i ,. the dissipation of fluctuations locally, as 

is related to the small-scale magnetic 

i,i2\ T 2 /l2\ i r • is usually the case in the presence of an 

energy as {ip ) = L£ H {b ), where L bH is J ^ 

• i i • -i i r i mi external forcing and dissipation. To obtain 

the typical horizontal scale of b. Thus, the 

,.rr . / n , this relation, we multiply the equation for ib' 

negative magnetic diffusion (second) term 

i i . , . . . , N , , , by ib' and then take the average 

and hyper-resistivity (third) term are closely 

linked through the small-scale magnetic 

energy <*>. Indeed, the negative diffnsiv.ty 2*<**> + ^<W>4W = + WW) 



and hyper-resistivity together conserve total 
(■?// 2 ), while shuffling the (^ /2 ) spectrum 
toward large scales. We now put equation 
(12) in the following form: 



Here, the integration by parts was used 
assuming that there are no boundary 
terms. We note that either when the 
stationarity condition is not satisfied or 
when boundary terms do not vanish, there 
, , 9 , „ , , . will be a correction to our results (Blackman 

K ' di{ilj) + L 2 bH diV 2 (ip) ' 1 1 & Field 2000). When (^ ,2 ) is stationary, the 

first term on the left hand side of equation 
(14) vanishes, simplifying the equation that 
where no summation over the index i occurs, relates (b 2 ) to Tj = (u^ 1 ) = €ij 3 {di(f)ip') to 
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the form: 



(b 2 ) = -[-TAW+B (ij'd, 



Note that in 2D MHD (B = 0), the flux is 
proportional to r](b 2 ). This balance reflects 
the conservation of (ip 2 ), which is damped 
only by Ohmic diffusion. The second term 
on the right hand side of equation (15) can 
be evaluated in a similar way as for Tj, i.e., 
by writing 

= (W-W), (16) 

and then by using equations (8)— (11) . 
Omitting the intermediate steps (see 
Appendix A for details), the final result is 

ty'M = rBoiUu 2 ) - Ub 2 )] ■ (17) 

Here 

& = J d 3 kk 2 0(k)/ J d 3 kk 2 H ^(k),(lS) 

6 = / d 3 kk^(k)/ y d 3 kk 2 H ^(k) ,(19) 

and k 2 H = k 2 + k 2 . If the characteristic 
horizontal and vertical scales of u are L v h 
and L vz , and if those of b are L^h and L^, 
then and £ 6 can be expressed in terms of 
these characteristic scales as: 



T 2 

Liz 



Cb 



T 2 

LI' 



Insertion of equation (17) into (15) gives us 

b 2 ) = ^[-TAm + r^B 2 (u 2 )]/(^l + ^B^l) 

Thus, from equations (13) and (21), we 
obtain 



r_ {u2) i + ^fe-^ + ^^l^r, . 



1 + 1 



\{Bn) 2 + i b B 2 - L Mj) 2 ] 



(20) 



where Jz = VxB ff and the integration 
by part is used to express <9j (ip) <9jV 2 (ip) = 
-(V 2 (^)) 2 = -(J) 2 < 0. Note the last 
term in the numerator and denominator 
in equation (22) comes from the hyper- 
resistivity. Equation (22) is the flux in 
3D RMHD, which generalizes the 2D 
MHD result (Cattaneo & Vainshtein 1991; 
Gruzinov & Diamond 1994). Several aspects 
of this result are of interest. First, in the 
limit as B -> and (B H ) -> ((J) -> 0), 
the flux reduces to the kinematic value 
Ti = —r)kdi(ip), with the kinematic turbulent 
diffusivity rjk = t(u 2 )/2. This corresponds to 
the 2D hydrodynamic result where the effect 
of the Lorentz force is neglected. The full 2D 
MHD result can be obtained by taking the 
limit B — > in equation (22), which will 
reproduce equation (1). This agrees with 
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the well-known result on the suppression of 
flux diffusion in 2D (Cattaneo & Vainshtein 
1991; Gruzinov & Diamond 1994). 

Another interesting case may be 
the limit (B#) — > 0. In fact, this limit 
can be shown to be consistent with the 
ordering of 3D RMHD as follows. First, 
note that 3D RMHD ordering (k z /kn ~ 
B H /B ~ e < 1) requires £ b B 2 ~ (B%). 
Since (B H ) 2 < (B%) ~ (b 2 ), we expect 
that i h Bl ~ (b 2 ) > (5^) 2 . Furthermore, 
L 2 bH (J) 2 ~ (L bH /L BH ) 2 (B H f < (B H ) 2 , 
where is the characteristic scale of 
(.Bff). Thus, the dominant term in the 
square brackets in the denominator of 
equation (22) is £ b B 2 ~ (b 2 ). That is, the 
effect of Bq seems to be stronger than that 
of (B H ) in 3D RMHD. 

Finally, to determine whether B 
enhances the flux or not, we note that 

— ^ in equation (22) can be taken to 
be zero, since the scales for b and u are 
likely to be comparable in this model, which 
employs unity magnetic Prandtl number. 
Then, we estimate the last term in the 



numerator, due to hyper-resistivity, to be 
r(b 2 )L 2 bH /(r]L 2 BH ) ~ (L bH /L BH ) 2 R m where 
£ b B 2 ~ (b 2 ) and (b 2 ) ~ (u 2 ) are used. If 
(L bH / L B h) 2 ~ -Rm 1 , this term will be of 
order unity. Note R m — ul/rj is the magnetic 
Reynolds number, with u and / being the 
characteristic amplitude and length scale 
of the velocity. Therefore, equation (22) 
indicates that the flux is reduced on account 
of the strong axial magnetic field B as well 
as the horizontal reconnecting field (B#). 
The above analyses will be used in §4.1 in 
order to estimate the effective dissipation 
and global reconnection rate. 

3. EDDY-DAMPED FLUID MODEL 

The analysis performed in the previous 
section introduced an arbitrary correlation 
time r that is assumed to be the same for 
both small-scale velocity and small-scale 
magnetic fields. Moreover, the quasi-linear 
closure is valid strictly only when the 
small-scale fields remain weaker than the 
large-scale fields. In order to compensate 
for these shortcomings, we now consider an 
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eddy-damped fluid model which is based a 
large viscosity (Kim 1999). In this model, the 
fluid motion is self-consistently generated by 
a forcing with a prescribed statistics as well 
as by the Lorentz force, without having to 
assume the presence of fully developed MHD 
turbulence, to invoke a quasi-linear closure, 
or to introduce an arbitrary correlation 
time for the fluctuating fields. This is the 
simplest model within which the nonlinear 
effect of the back-reaction can rigorously 
be treated. Even though this model has 
limited applicability to a system with a large 
viscosity, it could be quite relevant to small 
scale fields in Galaxy where v 3> rj. As shall 
be shown later, this model gives rise to an 
effective correlation time for the fluctuating 
magnetic fields that is given by the viscous 
time r u = l\ H jv, where l^u is the typical 
scale of the magnetic fluctuations in the 
horizontal plane (cf eqs. [22] and [32]). Thus, 
in comparison with the r approximation in 
the previous section, this model is equivalent 
to replacing r by r„ despite the fact that 
some of detailed results for the two models 
are not the same. 



3.1. Splitting of Velocity 

In a high viscosity limit with the fluid 
kinetic Reynolds number Re = ul/u < 1, 
the nonlinear advection term as well as 
inertial term in the momentum equation 
can be neglected. Then, the linearity of the 
remaining terms in the momentum equation 
enables us to split the velocity into two 
components; the first — random velocity - 
is solely governed by the random forcing, and 
the second — induced velocity — is governed 
by the Lorentz force only. Specifically, we 
express the total velocity u as u = v + v', 
where v and v' are the random and induced 
velocity, respectively, and introduce velocity 
potential O and 0/ as v = Vx0 o i and 
v' = Vx0/i. Then, the equations for these 
potentials are: 

= vV 2 <f) + F, (23) 
= z/V 2 0/ + B • VVV , (24) 

where the nonlinear advection term as 
well as the inertial term is neglected since 
Re < 1 is assumed. In equation (23), 
F is a prescribed forcing with known 
statistics. Instead of solving equation (23) 
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for 0o, we can equivalently prescribe the 
statistics of the random velocity O ( or v )- 
Therefore, we assume that the statistics of 
random component satisfies homogeneity 
and isotropy in the horizontal plane and 
homogeneity and reflectional symmetry in 
the z direction, respectively. Furthermore, 



3.2. Magnetic Field 

Both random and induced velocities 
are to be substituted in equation (2) to 
solve for the magnetic field. Notice that 
equation (2) then has a cubic nonlinearity, 
since the induced velocity is quadratic in B. 



we assume that it is delta correlated in time. We a S ain assume that the magnetic field in 

the horizontal plane consists of mean and 

fluctuating components, i.e., ip = (ip) + ip' 

5(ki - k 2 )5(ti - t 2 )0o^4i|h^2|^ e fluctuation is homogeneous and 

isotropic in the x-y plane and homogeneous 



The correlation function is then given by: 
(0o(ki,ti)0o(k 2 ,t 2 )) 



where (f>o(kiH,k lz ) is the power spectrum 
of 0o- Note that r o (0o) — J d 3 k<f)(k) and 
T o{v 2 ) — I d 3 kk 2 <i)(k), where r is the 
correlation time of v that is assumed to be 
short. 



and reflectionally symmetric in the z 
direction, satisfying the same correlation 
function as equation (10). 



To obtain equations for (ip) and (ip 2 ), we 
utilize the delta-correlation in time of v and 
On the other hand, the induced velocity iterate equation (2) for small time intervals 
can be constructed by solving equation (24) St. Specifically, we use (w i (ti) J B(t) j ) = 



for (pi in terms of B. This can easily be done f or t\ > t and 



v ~ 



0((5t) 



-1/21 



since 



in Fourier space as: 
0/(k) 



(vi(ti)Vj(t2)) oc -t 2 ) ~ 1/St, where 



vk 2 k 2 H 



St — t-i — t 2 . Then, for 5t 1, equation (2) 



where Bniik) = ieij^kjipik) is used. Note 
that the ip in the above equation contains 
both mean and fluctuating parts. 



^(t + St) 



t+st 



ip(t) + 5tr]V 2 ^(t) + / dt x [eijsdj^d^ih) + B c 
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\e lj3 f st dt ld t 2 {e lm3dl mdAdr r mdMt+ fe^fe^Si Twm&fi to 

the kinematic turbulent diffusivity and 



+ 



where ip and are to be evaluated at the 
same spatial position x. 

The mean field equation is obtained 
by substituting equation (26) in (27), by 
taking the average with the help of equations 
(10) and (25), and then by taking the limit 
St — > 0. The derivation is tedious and is 
outlined in Appendix B. Here, we give the 
final result 



4 x 1 2v 



1 



V 2 W 



Here To is the short correlation time of 
random velocity v and 



Vm 



— (v 2 ) — —G = rj k — —G , 
F 



\i = 



v 



r h 2 h 2 T 4 

F = J d 3 k^^(k) ~ ^f-G = 7 G , 

where r]k = tq(v 2 )/4 is the kinematic 
diffusivity; « = L 2 ^ and 7 = Lf H /L 2 bz = 
The above equation implies that the flux 
Ti = (liiV') is given by 

= -riMdiM+iidiV 2 ^). (29) 



backreaction. Note that the kinematic 

diffusivity i]k = r (f 2 )/4 now comes only 

from the random velocity, with r being 

its correlation time that can be prescribed. 

The backreaction term is proportional to 

(ip' 2 ), not (b 2 ) (cf. eq. [11]) and inversely 

proportional to the viscosity v. It is because 

the cutoff scale of the magnetic field l v 

is smaller than that of the velocity l v in 

tm^n^ile^so that for a larger v, there are 
v 

magnetic modes over a larger interval of 
(28) 

scale I between l v and l v (i.e. l v < I < l v ) 
where the velocity is absent due to viscous 
damping. That is, the induced velocity 
(Lorentz force) cannot be generated on this 
scale (l v < I < l u ) due to viscous damping, 
thereby weakening the overall effect of 
backreaction (see eq. [48]). Now, the last 
term in equation (29) is the contribution 
from the hyper-resistivity fi. It is interesting 
to see that ji is inversely proportional to 
L\ z and thus vanishes as L hz — > 00 (or 
7^0) which corresponds to the 2D limit. 
Therefore, in this eddy-damped fluid model, 
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the hyper- resistivity term vanishes in two In a stationary case, equations (28), (30), 

dimensions. It should be contrasted to the and (31) lead us to the following expression 

case considered in the previous section where for the flux: 

the hyper-resistivity, being proportional _ 1 + ^-B' 2 {^ b — £„) + -^£ v Bq 



= —t( v ) : z r . ^ , — TiZ-bm: 



(V> /2 >, survives in 2D MHD limit (see eq. 4 1 + £ [&B 2 + \{B H ) 2 - 7 (J) 2 ] 

[12]). where Jz = VxB^, and d& 2 (ip)di(ip) = 

_(V 2 (^)) 2 = -(J) 2 < is used. When the 
characteristic scales of fluctuating velocity 
and magnetic field are comparable, or when 



For use later, we solve equation (29) for 
(b 2 ) yielding 



2 = Tj + ^dM 
\ ' _|_ «ig.y2^ ' ^ ' only the ratios of vertical to horizontal scales 

where again the summation over the index i of the fluctuating velocity and magnetic 
is not implied fields are comparable, can be taken to be 

equal to simplifying the above expression. 

3.3. Stationary Case: d t (^' 2 ) = 

It is worth considering a few interesting 

The additional relation between the flux ,. ., ( / QO \ t?- , • , u r ., 

limits of equation (3zJ. rirst, m the limit 

Tj and magnetic energy (b 2 ) is obtained for D n , D n ,. / QO x 
* to 6j \ / ^ Q o and i% — >■ 0, equation (32) again 

the case of stationary (^' 2 ). To derive an ^ 2D hydrodynamic result with 

equation for (^ 2 >, we multiply equation (27) thfi kinematic diffusivity % = To <^)/ 4 . 

by itself, take average, and then take the Thg ^ ^ ^ Q le&ds tQ m MRD ^ 

limit of 5t — > 0. After considerable algebra , . r +u + u i + 

to where the suppression of the turbulent 

(see Appendix B), we obtain the following diffusion ^ from ^ 3D RMRD 

equation ^ e dominant suppression in the flux comes 

d t (^) + d t {i>) 2 - 2 V [-((d^f) + (^)V 2 (^)] f ^ m ^o $^)& =$ljs discussed in §2.3. 

We note that the last term in the 



where 



2 2 numerator and denominator is due to 



G = [ d 3 k^ij(k) ~ i^G = , . ... 

7 k 2 L z hH the hyper-resistivity, which comes with a 
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multiplicative factor 7 = L\ H jL\ z <C 1. 
Therefore, the effect of hyper-resistivity 
can be neglected as compared to other 
terms in equation (32). Since 7 — > in 
2D MHD, there is no contribution from 
the hyper-resistivity to the flux in 2D in 
this model. The estimate of the effective 
dissipation in this model is provided in §4.2. 

It is very interesting to compare 
equation (32) with (22). We recall that 
in order to derive equation (22), the same 
correlation time r was assumed for both 
fluctuating magnetic field and velocity, which 
appears in front of the mean magnetic fields 
B and (ip) in equation (22). In contrast, r 
in equation (32) is the correlation time of 
the random component of the velocity, which 
can be arbitrarily prescribed. Moreover, r 
in front of mean magnetic fields in equation 
(22) is now replaced by viscous time scale 
t u = k/v = L\ H jv in equation (32). The 
latter represents the viscous time scale across 
the typical horizontal scale of fluctuating 
magnetic fields. Thus, as noted at the 
beginning of this section, this viscous time r u 
replaces r in the quasi-linear closure, which 



was assumed to be a parameter. 

4. RECONNECTION RATE 

In previous sections, the flux Tj was 
derived by using a quasi-linear closure and 
an eddy-damped fluid model. Assuming the 
flux Yi has a form proportional to di(ip) 
in both cases (see eqs. [22] and [32]), it 
can be expressed in terms of the effective 
dissipation rate (or, turbulent diffusivity) 
r] e ff as follows: 

= -VeffdiOj). (33) 

Upon using equation (33), the mean field 
equation (4) then becomes 

dM) = (^ + r ?e// )V 2 W=r ?T V 2 (^X34) 

where r] T = r\ + r) e ff is the total dissipation 
rate of the mean field. The effective 
dissipation rate is the quantity that 
represents the overall decay rate of a large- 
scale magnetic field due to both small-scale 
motions and magnetic fluctuations. That 
is, the dynamical system consisting of 
both small and large scale fields can be 
represented by the evolution of a large-scale 
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field only when the effect of small-scale fields 
is absorbed in this turbulent coefficient. 

In order to determine a global 
reconnection rate, we now invoke the original 
SP type balance equations and use the 
total dissipation rate in place of the Ohmic 
diffusivity (see §2): 



4.1. 



V r = 



VA 



(35) 



Note that we have neglected a multiplicative 
correction factor to the reconnection 
rate in the eddy-damped model since its 
dependence on v is weak with 1/4 power 
(for instance, see, Biskamp 1993). In the 
following subsections, we assume £ v = 
for simplicity and estimate the reconnection 
rate via equation (35). Then, we briefly 
comment on the implication for reconnection 
assuming 'Alfvenic turbulence', as Lazarian 
and Vishniac (1999) did. 



Using the Quasi-linear Result 



The effective dissipation rate follows 
from equations (22) and (32): 



Veff 



1 _|_ TL hH t R 2| 9 t V 2 W 



1 + T - 

V 



after using £ v = £ 6 . As shown in §2.3, the 
dominant term in the square brackets in the 
denominator of equation (36) is £ b B 2 ~ (b 2 ), 
and the second term in the numerator is of 
order unity if L 2 BH jL 2 hH ~ R m . In that case, 
r] e ff is roughly given by 

1 1 



Veff 



Vk- 



Vk- 



l + r<6 2 )/77 l *l + 2R m (V) 
where r] k = r(u 2 )/2 is the kinematic value of 
turbulent diffusivity in 2D and R m = rjk/rj. 
In contrast to the 2D MHD result (eq. [1])), 
the equation (37) reveals that the effective 
diffusivity in 3D RMHD is more severely 
reduced as {b 2 } > {B H f (= (B) 2 ). To 
determine the leading order contribution in 
equation (37), we need to estimate (b 2 ). To 
do so, we substitute equations (33) and (37) 
in (13) and use < Lbh to obtain: 



T 

where R m = rjk/r) = t{u 2 )/2t] is used. 
We note that (b 2 ) > is guaranteed since 



1 - 



2R r , 



(38) 
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(b 2 ) > (B H ) 2 (implying R m > 1) was 
assumed to derive the above equation. Thus, 

r(b 2 ) 



V 



That is, for R m » 1, T(b 2 )/i] > I. Insertion 
of the above equation in (37) then gives us 

1 rj 



Veff ~ Vk 



2R m 2 



(39) 



In other words, to leading order, the effective 
dissipation rate is just that given by Ohmic 
diffusivity! Therefore, by inserting equation 
(39) into (35) with t]t — rj + rj e ff, the 
reconnection rate is found to have the 
original SP scaling with rj, i.e. 



V v aL/i] 



(40) 



It is interesting to contrast this result 
to the 2D case where B = 0. In that 
case, the dominant term in equation (36) 
is (B H ) 2 , with r} eff ~ i] k (u 2 ) / R m (B H ) 2 ~ 
r)(u 2 ) I (B H ) 2 ~ 7]u 2 /v 2 A > i], where u is the 
typical velocity. Therefore, in 2D, the global 
reconnection rate becomes 



V A u 



(41) 



which is larger than SP by a factor of 
magnetic Mach number M A = u/va- Note 



that the reduction in the effective dissipation 
of a large-scale magnetic field is more severe 
in 3D RMHD than in 2D MHD by a factor 
of (u 2 )/(B H ) 2 ~ {u 2 )/v\. 



4.2. Using the Eddy-Damped Fluid Model 

Result 

For an eddy-damped fluid model, 
equation (32) yields: 



Veff 



T 0/ r 2\ " ' ^"~ U *W (AO 



1 + 2 ^B 2 



after assuming £ v = £ b . We recall that the 
contribution from the hyper-resistivity comes 
with a multiplicative factor 7 = L 2 H /L 2 Z 1 
(vanishing in the 2D MHD limit) and thus 
can be neglected as compared to other terms 
in equation (42). Then, a similar estimation 
as in §4.1 simplifies equation (42) to 



Veff 



(43) 



where rjk = r (f 2 )/4 is the kinematic value of 
the turbulent diffusivity in 2D and k = L 2 H . 
To obtain the leading order behavior of 
equation (43), we estimate (b 2 ) with the help 
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of equation (30) to be is (B H ) . Thus, r/ e ff ~ r] k (v 2 ) / R e R m (B H ) ~ 

„ v r](v 2 )/R e (B H ) 2 ~ r]u 2 /R e v 2 A > rj, where w is 

(fc 2 > ~ -^-(2i? m -l), (44) 

K the typical velocity. Therefore, in 2D, the 

where i? m = By inserting equation global reconnection rate becomes 

(44) in (43), we obtain 1 va u 

v r ~ —=—= — , (49) 

VR e Jv a L/t] v a 

Veff ~ 2^ ~ 2 ' where u/va = M A is the magnetic Mach 



™ ,. , number. In comparison with equation (41), 

Thus, the reconnection rate is again given by v ' 

the global reconnection rate in this model is 

^ ~ ^ VaL/7] ' ^ tllUS lalger 111 2D Umit ( reCaU R e <1 )- 

i.e., SP scaling with 77 persists! 

It is interesting to estimate (b 2 ) in 4 - 3 - Alfvenic Turbulence 

equation (44) by using 

In Alfvenic turbulence (Goldreich & 
T- = ( w 2 ) 1 - ^ (t; 2 ) 1 (4^ridhar 1994; 1995; 1997), the equipartition 

K \fW)L b H \/W)L bH R m Re 

between (b ) and (u ) is assumed from the 

where R e y/Jv )L bH /v is the fluid start. It is to be contrasted to the present 

Reynolds number. Thus, equation (44) , . . , . , , , . . 

analysis m which the relation between 

becomes ^ an( j j e ^ e q Ua tions (38) and (49), 

(b 2 ) ~ (v 2 )— (2 —J (48) f°ll° ws from the condition of stationarity 

R e \ Rm J 

of (ip' 2 ) in the presence of B and (B H ). 

The above equation clearly demonstrates . , . . /ooN . , 

As can be seen from equation (38), m the 

that (b 2 ) > (v 2 ) for our model (R e < 1) when ... ... 

quasi-lmear closure with unity magnetic 

R m > 1, as pointed out near the end of §3.2. .... 

rrandtl number, exact equipartition is 

Finally, we note that in 2D limit with possible only for r] = 0. In the eddy-damped 

Bq — > 0, the dominant term in the square fluid model, exact equipartition can never be 

brackets in the denominator of equation (42) satisfied since the assumption R e < 1 implies 
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(b 2 ) > (v 2 ) when R m > 1 (see eq. [48])! 
Therefore, in general, stationarity of {i/j' 2 ) 
and exact Alfvenic equipartition cannot be 
simultaneously achieved. In other words, if 
Alfvenic turbulence is assumed, (vp' 2 ) cannot 
be stationary; if (ip' 2 ) is stationary, the 
turbulence cannot be in a state of Alfvenic 
equipartition. 

We easily confirm this in 2D MHD by 
quasi-linear closure. The exact equipartition 
({u 2 — b 2 ) = 0) implies that the flux in 
equation (12) is given by hyper-resistivity 
only: = -r(f 2 )d l V 2 (^)/2. Then, 
if we were to impose the stationarity 
of (ip' 2 ), equation (15) would indicate 
(*fj' 2 )Td(J)(B H ) = V (b 2 ). Thus, 

where Ib and If, are the characteristic 
scales of (B#) and b, respectively. Since 
(B H ) 2 /(u 2 ) ~ l/R m (with (b 2 ) ~ (u 2 )) and 
(h/lb) 2 ~ l/Rm in 2D MHD, the relation 
(49) (for stationarity) cannot be satisfied. 



5. CONCLUSION AND 
DISCUSSIONS 

In view of the ubiquity of turbulence 
in space and astrophysical plasmas, 
magnetic reconnection will likely occur in 
an environments with turbulence. On the 
other hand, the reconnection itself generates 
small-scale fluctuation, feeding back the 
turbulence. Thus, it is important to treat 
these two processes consistently, accounting 
for the back reaction. Although LV argued 
that the local reconnection rate can be 
fast, they basically neglected the dynamic 
coupling between small and large scale fields, 
therefore leaving the issue of the global 
reconnection rate unresolved. The coupling 
between global and local reconnection rates 
should be treated self consistently. The aim 
of the present work was to shed some light 
on this issue by taking the simplest approach 
that is analytically tractable. 

Our main strategy was to self- 
consistently compute the effective dissipation 
rate of a large-scale magnetic field within the 
current sheet by using stationarity of (ip' 2 ) 
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and then use the effective dissipation rate in 
SP type balance relations to obtain the global 
reconnection rate. To avoid the null point 
problem associated with a 2D slab model, 
we considered 3D RMHD, within which 
we can solidly justify the incompressibility 
of the fluid in the horizontal plane. To 
facilitate analysis, two models (methods) 
were employed, one being a quasi-linear 
closure with r approximation and the other 
eddy-damped fluid model. 

The effective dissipation rate r) e ff that 
we obtained generalizes the 2D MHD result 
(Cattaneo & Vainshtein 1991; Gruzinov & 
Diamond 1994). The quasi-linear closure 
predicted i] eff ~ rj k /(l + 2R m (b 2 ) / (u 2 )) ~ 
?7/2 (see eqs. [37]-[39]). A similar result was 
obtained in the eddy-damped fluid model 
with ri eff ~ 7fo/(l + R m R e (b 2 )/(u 2 }) ~ 77/2 
(see eqs. [43]-[45] and [47]). 

The 2D result can simply be recovered 
from our results on the flux by taking 
the limit B — > 0. In that limit, 
Veff ~ W(! + Rm(B H ) 2 /(u 2 )) according 
to the quasi-linear closure, consistent with 



previous work. In the eddy-damped fluid 
model, r] e ff ~ r) k /(l + R m R e (B H ) 2 / (u 2 )). 

Since the effective dissipation rate 
i] e ff was found to be the same in both 
models (in 3D RMHD), the global 
reconnection, obtained by invoking SP 
balance relations, was also the same with the 
value v r ~ va/ ^/vaL/v in both models. This 
result indicates that the global reconnection 
rate is suppressed for large R rn as an inverse 
power of Rm 2 such that the original SP 
scaling with r\ persists. Again, this persistent 
77 scaling results from the reduction in the 
effective dissipation rate of a large-scale 
magnetic field for large R m mainly due to a 
strong axial magnetic field, with the effective 
dissipation rate r) e ff ~ rj. 

Furthermore, in the 2D limit, the 
quasi-linear closure yielded the global 
reconnection rate v r ~ (va/ \/vaL/v)(u/va), 
which is enhanced over SP by a factor of 
Ma — u/v A (note that M A can be large). In 
contrast, the eddy-damped fluid model gave 
v r ~ \fR~ e 1 {y A / \fv A L/ri){u/v A )- 

The implication of these results for 
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the LV scenario is that no matter how fast 
local reconnection events proceed, there is 
not enough energy transfer from large-scale 
to small-scale magnetic fields to allow 
fast global reconnection. Therefore, global 
reconnection cannot be given by a simple 
sum of the local reconnection events as 
LV suggested. We emphasize again that 
the {ip' 2 ) balance played the crucial role in 
determining the global reconnection rate 
consistently. Alternatively, an accurate 
calculation of the global reconnection 
rates requires that (global) topological 
conservation laws be enforced. 

The reduction in the effective dissipation 
in 2D is closely linked to the conservation 
of mean square magnetic potential. In 
3D RMHD, the mean square of parallel 
component of potential is no longer an 
ideal invariant due to the propagation of 
Alfven waves along a strong axial magnetic 
field. Nevertheless, the conservation of mean 
magnetic potential is broken only linearly, 
which turned out to introduce additional 
suppression factors, as compared to 2D. The 
interesting question is then how relevant 



these results would be in 3D. The mean 
square potential is not an invariant of 3D 
MHD. However, its conservation is broken 
nonlinearly, unlike 3D RMHD. Therefore, 
the effective dissipation in 3D MHD may 
be very different from that in 3D RMHD, 
with the possibility that the former may not 
be reduced, at least, in the weak magnetic 
field limit (Gruzinov & Diamond 1994; 
Kim 1999). Moreover, in 3D, there is a 
possibility of a dynamo, which brings in 
an additional transport coefficient (the a 
effect) into the problem. Some insights into 
the problem of effective dissipation of a 
large-scale field in the presence of a dynamo 
process might be obtained by considering a 
simple extension of the present 3D RMHD 
model by allowing a large-scale dynamo 
in the horizontal plane. Recall that this 
possibility was ruled out in the present 
paper by assuming isotropy in the horizontal 
plane and reflectional symmetry in the axial 
direction, with no helicity term (i.e., no 
correlation between horizontal and vertical 
component of fluctuations). 

Considering some of limitations of 
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the two models that were analyzed in 
the paper, such as the r approximation, 
quasi-linear closure, low kinetic Reynolds 
number limit, etc, it will be very interesting 
to investigate our predictions via numerical 
computation. The stationarity of {ip' 2 ) 
can be maintained as long as there is an 
energy source in the system, such as an 
external forcing. By incorporating the 
proper ordering required for 3D RMHD, one 
can measure the decay rate of (B#) to check 
our predictions for r) e ff ~ r) (see eqs. [40] 
and [46]). Ultimately, a numerical simulation 
with a simple reconnection configuration 
should be performed to measure a global 
reconnection rate as a function of R m as well 
as Bq and {B H }. It will also be interesting 
to investigate non-stationary states such as 
plasmoid formation (Forbes & Priest 1983; 
Priest 1984; Matthaeus & Lamkin 1986). 

We thank E. Zweibel for bringing this 
problem to our attention and for many 
interesting discussions. We also thank E.T. 
Vishniac and A.S. Ware for stimulating 
conversations. This research was supported 



by U.S. DOE FG03-88ER 53275. PH. 
Diamond also acknowledges partial support 
from the National Science Foundation 
under Grant No. PHY99-07949 to the 
Institute for Theoretical Physics at U.C.S.B., 
where part of this work was performed. 
E. Kim acknowledges partial support from 
HAO/NCAR where part of this work was 
completed. 

Appendix A 

In this appendix, we provide some 
of steps leading to equations (12) and 
(17). First, to derive equation (12), we let 
r, = rj 1} - if \ where if > = e lj3 (<W> and 
if ^ = 6ij3((f)diil>') , and begin with rf . 

= e ij3 J rf 3 A; 1 rf 3 A; 2 iA; lj (0(k 1 )^ / (k2))exp{i(k 1 + k 2 )(- 

After inserting equation (8) in (Al) and 
using equation (11), we can easily obtain 

r t (1) = -ire ij3 ei m3 j ' rf 3 A; 1 rf 3 A;A; li A; im A; i 0(ki)(V'(k))e ikx + t 
= —9^)5,1 J d^iA^ki) = — (u 2 )^) . 

where (u 2 ) = j d?kikl<f)(lLi). To obtain 
the last line in equation (A2), we use the 
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following relations A bit of algebra then gives us 

d 3 kk j k n ^(k) = U jm J d 3 kk 2 0(k) , rf) = L [-(b^dii^x)) - (^ 2 ><9,V 2 (?K^.]7) 
d 3 kk j k z 4(k) = 0, (A.3) Thu§) from equationg (A3) and (A7); we 

which follows from the isotropy of in the obtain equation (12) in the main text. 
x-y plane, and reflectional symmetry in the g econd; tQ deriye equ&tion ^ wg 

z direction. . , , ,. f 

again compute the correlation function 

The second part, rf\ is calculated in a on the right hand side of equation (16) 
similar way. in Fourier space. The first term can be 

rewritten as: 

= e ij3 J d 3 k 1 d 3 k 2 ik 2j (8(f)(k 1 )^(k 2 )) exp {^A^Jk^j} d 3 k 1 d 3 k 2 ik lz ((j)(k 1 )54>'(k 2 )) exp {i(ki + k 2 ) 

We insert equation (9) in (A4) and use (10) Then, inserting equation (8) in (A8) and 
to obtain using equation (11) give us 

-iB J tfhkuhrfik!) (Srj/d^) = r[Jd 3 k 1 k lz k lz B ^(k 1 )-e lm3 Jd 3 k 1 d 3 kk lz k 

(k + k2)2 [kmhikj + k 2m hk^ ^# (y 2 )#^||(k 1 ) = tBoUu 2 ) , 

Since (ib) has a scale much larger than ib', i ^ ■ i ^ iia\ 

x ' where the isotropy and equation (18) were 

k 2 k in the second integral on the right , , , , . , , n . , 

used to obtain the last line. Similarly, the 

hand side. We thus expand the integrand , , , , . , , . , , 

second term on the right side of equation 

of this second term and use the following ■ -i i , . , /• • \ i 

(16) is easily calculated (m rouner space) by 

isotropy relations: . ™ . 

using the isotropy condition. The result is 

/ d 3 kkikjkik m il)(k) = - (5^5^ + 5 u 5 jm + 5 im Sji) / d 3 kk 4 ijj(k) , 

J 8 Thusf equations (16), (A9), and (A10) yield 

d 3 kkikMk) = 0. ,. /^(A.61 . , , 

J rv ' equation (17), v m the mam text. 



r (2) ire 
1 i — lT€ ij3 



+Qm3 / d k 2 d ke 



-27- 



Appendix B 

In this Appendix, we provide some of 
intermediate steps used to obtain equations 
(28) and (31). For the mean field equation 
(28), we first take the average of equation 
(27) 



{vfatjvjfah)) = 5{h - t 2 )%T L (r = 0) 
so that T L (0) = t (v 2 )/2 = 2rj k . Here 
Tq is the short correlation time of v 
and i]k = r (f 2 )/4 is the kinematic 
diffusivity. (vi(x)^(x)) is obviously re- 
lated to 0o by (d i (f) (x,t 1 )d l <f) (y,t2)) = 



Sa(vj(x.,t 1 )v j (y,t 2 )) + (v i (x.,t 1 )v l (y,t 2 )). By 
(tp(t + 8t)) - (ip(t)} - 5trjV 2 (i/j(t)) = h + Vj = -€i j3 di<j>o and (0 o (ti)^(t)) = 0, 

Is is determined to be: 



where 



ft+ st h = -5tT L (0)V 2 (^). 

h = dti hjsdjipi^di^iih)} ~ Ste^diidjipit^iit)) = 5tdi\ , 

1 t+ st h represents the kinematic turbulent 

h = 



(B.4) 



dt lBo d z (Mh)) - stB d z (Mt)) )diffusivity . Nextj to compute /2> we take the 

i rt+St 

+B a i Mti)d j ,Mt2)) , and then take the tW e ' Upon neglecting 

d z {ip) ~ 0, one can easily show that I 2 = 0. 
where A, = e ijZ (d^{t)h{t)) and the smooth Finally? h contains the backreaction as well 
variation of the induced velocity 0/ in time ag hyper-resistivity. To evaluate this term, 
was used to approximate the time integrals 
in Ii and I 2 . To compute the averages, 
it is convenient to express the correlation 
function (24) in terms of v in real space as: 



we insert equation (26) in Aj to obtain 

A, = eijsidtfWMt)) 



v 



j3;„/ ik'-x 



— Uj3^im3 / d k 2 d k e 



(k + k') 2 (k^ + k^ 



(ui(x,*i)uj(y,i 2 )) = S(t 1 -t 2 ) 



here dTj 



T L (r H ,r z )S i:j + r H 



where r = y — x and r# is the hor- 
izontal component. Note that the 
above relation implies that at r = 0, 



dr H 

Pjlm 



r H ir Hl 



r 



H 



3) 



-kj \k m k l k'^ l + /c;/c m /c^-] . 



For notational convenience, we introduce 
q = k# so that g 3 = 0. Since the 
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characteristic scale of (ip) is much larger where 
than that of ip', k! <C k in equation (B5). 
Thus, we expand the integrand of equation 



than that of ip', k! <C k in equation (B5). rt+st 

rt+St 



dt 1 dt 2 {ip{t) [di4>o{ti)ei m 3dj [d m ip(t)di<f> (t2) 

the isotropy and homogeneity of tp' in the f t+5t 

' h = ^(^(f)[6 i ^(f)^(i 1 )+Bo^ / (t 1 )]>E l 

x — y plane. The latter implies equation Jt 

(A5) (recall q = k H ) and also the following where J 31 = €ij 3 {ip(t)djip(t)di(f)j(t)) and 
relations J 32 = B (ijj(t)d z (f)i(t)). 

J d 3 k m q r k n = J d 3 k m q r q n , First ' J i can easil y be computed by 

/I r using the correlation functions as 

d 3 k m k z k z = -5 jX J d 3 kq 2 k 2 .(B.6) 

Ji = St 

Then, a fair amount of algebra reduces 



T L (0) [<6 2 > + (B H ) 2 ] -Blj d 3 k 2 z <m0) 



equation (B5) to Next, J 2 can be computed upon substituting 

1 f k 2 1 Wuati$ru426) and then splitting average by 

a, = ~m / d^m - -9,v 2 ^) 

^ J k " fusing (^(i)^(tx)) = 0, with the result 

= -^^)-^V 2 (V). (B.7) 

J 2 = 5tT L (0)[-(6 2 ) + WV 2 W}B.ll) 

Note that there is no contribution from the 

first order term. By inserting equation (B7) For J 3, one can first show J 31 = due to 

into (Bl), by dividing both sides by St, and ^otropy To compute J 32 , we substitute 

then by taking the limit of St -> 0, we obtain equation (26) and use (0,) = to obtain 
equation (28). 



Next, to derive equation (31), we 

x(^'(ki) 

multiply equation (27) by ip and then take 
average to obtain the following equation: 



/k 
d 3 hd 3 kexp {i(k x + k 2 ) -xj-p^ 

ok z k 2 H iP'(k) + ie ij3 J d 3 k'*P(k - k')(k 



k 2 

mt + 5t } ) - ( ^ 2 (t)) - 2nstw)v 2 m) = tmJ Mkl exp Wkl + k2) ' x} 5fepi Q ^ (ki: 
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where Q tj = —kij(k + ki)i(k H + k m ) 2 - 
ku(k + ki)jkf H . By using the definition 
of G (see immediately after eq. [31]) and 
tijzQij = -tijsitfj + 2k H ■ k H )kikij, equation 
(B12) becomes 



J32 — 



.El 
v 

Bo 
v 



B G - ie ij3 I d 3 k'dV-*-!p- [k% + 2k x {k' - k)i] h{k' - k)^{-k + k')(V(k')> 



B G - e m d j( i 6 k' / d 6 ke 



j3uJk'-xk z ki 



k 2 



-1 + 



2kik\ 



^(-k + k')^(k')) 



(B.13) 



Now, since k' <C k, we expand the integrand 
of equation (B13) to second order in k'/k, in 
order to show that there is no contribution 
from the second term in equation (B13) 
(to this order). Therefore, J 32 = —BqG/u. 
Inserting J 1; J 2 , and J3 in equation (B8), 
dividing by St, and then taking the limit 
St — > finally yields equation (31). 
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Fig. 1. — Sweet-Parker 2D slab configura- 
tion. A and L are the thickness and length of 
the current sheet; ±B are reconnecting mag- 
netic fields; v r and vq are inflow (reconnec- 
tion) and outflow velocities. 

Fig. 2.— Configuration in 3D RMHD. B is 
a strong axial magnetic field pointing in the z 
direction, and ±B# are reconnecting (large- 
scale) magnetic fields in the x-y plane. Panel 
(a) shows the projection in the x-y plane, and 
Panel (b) in the y-z plane. 




(b) 



